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I. Introduction 

The role of rotation and magnetic fields in cosmical fluid dynamics 
is well known. Recently, some attention has been given by several 
authors [1-7] on the interaction of Coriolis forces and Lorentz forces 
on flow phenomena. 

It is generally accepted [l, 2 ] that hydromagnetic flow in the 
earth’s liquid core is somehow responsible for the main geomagnetic 
field, and therefore a theory of the dynamics of core motions is 
required in order to understand a plausible theory of the earth T s magnetic 
field. It has been suggested that the near coincidence between the 
geographic and geomagnetic poles is the result of the strong influence 
of Coriolis forces, due to the earth’s rotation, on motions in the core. 

So much of meteorology depends ultimately upon the dynamics of 
a revolving fluid. The large-scale and moderate motions of the atmosphere 
are greatly affected by the rotation of the earth. Kelvin is said 
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to have pointed out [8] that rotation confers on a fluid certain 
properties resembling those of an elastic solid, and a rotating 
fluid can transmit "waves . In the case of an infinite liquid, rotating 
as a rigid body about an axis, the amount of energy possessed by the 
liquid is infinite and It is of great interest to know how small 
disturbances propagate in such a liquid. 

To understand some of these phenomena, it will be interesting 
to study the flow of a rotating fluid around elementary bodies. 

We consider here small oscillations of a sphere in a compressible, 
viscous, electrically conducting and rotating medium in the presence 
of a uniform magnetic field. In addition to the above mentioned 
interests, this problem may have sane applications in connection with 
the transmission of sound by fog [see 9; P- 6591* 

The classical problem of the oscillation of a sphere in a viscous 
fluid was first considered by Stokes [ 9 , P* 643]* The drag force 
experienced by the sphere is (in dimensionless form)* 


X 

o 


-6 tt u 

s 



„ a Uc 

TT b 

d t 




(1) 


It is the purpose of the present investigation to determine the 
effects of the magnetic field, compressibility, coriolis forces and their 


*The dimensionless quantities and parameters in (l) are explained 
in section II. 
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mutual interaction on the drag and other physical variables. Because 
of the motion of the fluid in the magnetic field, an associated electrical 
field is produced which, according to Ohm's law sets up electrical currents 
in the fluid if the latter is a conductor. The interaction of these 
currents with the magnetic field then produces a body force which must 
be included in the Navier- Stokes equations for the motion of the fluid. 

The effect of this body force is to inhibit the motion of the fluid 
across the lines of force. The viscous effect gives rise to viscous 
dissipative waves, the magnetic effect is responsible for Alfven waves 
while the compressibility and rotation produce sound and the so-called 
Taylor waves respectively. Their mutual interaction will make the 
situation even more complicated. In the next section, we make several 
assumptions to make the problem mathematically more amenable. 



II. Fundamental Equations and Formulation of the Problem 


Let the sphere oscillate along the axis of the rotating fluid 
which is taken to he the x-axis. The origin is at the mean position of 
the center of the sphere. Let the uniform magnetic field he H°i. 

The equations of motion of a compressible, viscous, unbounded fluid 
rotating with a constant angular velocity Qi, referred to a rotating 
frame of reference are (in m.k.s. units) 


d v , 

P - + 2pQ)(v + pfl)((nXr) = - Vp+^ curl H X H +- 


( 2 ) 


+ p v V2 y + 21 V div v , 


where p is the density, v the velocity of the fluid, H is total magnetic 
field, p is the pressure, r is the radial coordinate: r 2 = y 2 + z 2 , 
v the kinematic viscosity (assumed constant) and p the magnetic 
permeability (assumed constant). 

The equation of continuity is 

div (p v) + = 0 . (3) 


The physical equation (the conduction of heat is neglected) is 
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(where 7 denotes the ratio of two specific heats) which can he 
written as [ 9 * p. 65^1 

p = p + p c^ s, (it) 

where c is the velocity of sound in the absence of viscosity, p°, p° 
represent pressure and density in the undisturbed state and s denotes 
the condensation 

P = p° (1 + s) . (5) 

Maxwell's equations with the usual notation for electromagnetic 
quantities are 

(i) curl H = J, (ii) div H = 0, 

SI (6) 

(iii) curl E = - p -yj- , (iv) div E = 0, 

(v) J = cr [e + p v X H] , 

where the electrical conductivity o is treated as constant. 

From (i), (iii) and (v) we obtain 

|| - curl (y X H) = 52 5 • 


( 7 ) 
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As in the classical problem, the convective terms in (2) will henceforth 
be neglected. This is justified if 

R = a |v |/v « 1 , (8) 

s s 

where v is the velocity of the sphere and a is its radius, 
s 

Further linearization of (2) and (3) is possible by means of (5) 
which yields 

p° -g-* + 2 p° A ^ X r)f = - ^ p + p curl H X H ... 

; " »> 
+ p° v V 2 v -f — — ^ V div v 
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and 


div v + = 0 . (10) 

These equations will be now reduced to a dime ns ionless form by referring 
the length to a, the radius of the sphere, the velocity to oX (\/2tt 
is the frequency of the oscillation), the time to l/X, the magnetic 
field to H°, the electric field to jj.aXH° , the current density 
J to H°/a and the pressure to p° a 2 X 2 . The equations (6), (7) then 


become 
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(i) curl H = J, (ii) div H = 0, 

(iii) curl E = - -|-S, (iv) div E = 0, 

(v) curl H = R [E + v Y H] 

— m — — n — 

and 

v2 H = R m [_ |-£ - curl (v X H) j 

where 

(Magnetic Reynolds number) = a 2 1 p cr 
In addition, we introduce the following parameters: 

R g (Reynolds number)* = a 2 X/v 


*It should be noted that we have used aX as velocity in the 
definition of the Reynolds number. There still is another important 
dimensionless parameter cr/a, where ot is the amplitude of the 
oscillatory motion of the sphere. It is assumed that this quantity 
is small compared to unity. 

The velocity of sphere is v s i = a l e^^ 3., or in dimensionless 
form, u s i. = u/a e^ _i. 

The Reynolds number for the flow might also have been defined 
as 


R = a <y X/v (=R a/'a). 

s ^ 


It is necessary that 


R « 1 

s 


y 


(11) 

( 12 ) 
(13) 

(14a) 


whereas the requirement for Re is less severe. 
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3 (magnetic pressure number) 


_ M- 


H 


o 2 




\2 


(l4b) 


the ratio of the magnetic pressure to dynamic pressure, and 


6=2 o/x 


(14c) 


the inverse of the so-called Rossby number R q = X/£Q) = l/6 . In terms 
of these parameters, equations (9) and (10) reduce to 


5 I 
s t 


+ 6 i X v 


3 curl H X H + v2 v-(c 2 + v— ^ v s 

± ' e ? e 


(15) 


and 


S s 

div v = - -5- 


( 16 ) 


where velocity of sound c is now in non-dimensional form referred to 
, and if we assume that 

R » 1 , (17) 

o 

so as to neglect second and higher powers of 6, the centrifugal force 
in (9) being of 0 (6 2 ) is dropped* in (15). 


^Notice that for large values of radial distance, our neglecting 
centrifugal force may not he justified even though S 2 is negligible. 

In the Incompressible case, it presents no difficulty since it can be 
absorbed with pressure term. However, in Compressible case, this cannot 
be done. But since we are primarily interested in the effects of 
Coriolis forces on the motion, following the standard procedure in such 
cases (see for example Hide and Roberts i960 [2] and Talwar 1964 C 5 3 ) 
centrifugal force may be ignored. 


- 9 - 


The boundary conditions on v are v = i on the surface of the 
sphere and is zero at infinity. If the conductivity o 1 of the sphere* 
is comparable to the conductivity o of the fluid, we require that 
H should be continuous across the surface of the sphere, while at infinity, 
H = i.** 


We next make two important assumptions: 

(1) The magnetic Reynolds number R^ is a small parameter. 

We can then find a perturbation solution for small R . This is a 

m 

valid assumption in most practical problems and this technique has been 
adopted in several discussions (see for example, Ludford [lO] and 
Tamada [ll]). 

(2) The magnetic pressure is comparable to dynamic pressure, 
i.e. we assume 3 = 0(l). 

We assume an expansion in powers of R^: 


v = v + R Vi + 

— — o m — 


H = H + R Hx + 

— — o m — 

E = E + R Ei + 

— — o m — 


(18) 


S = S + R Sn + 

o m x 


*The discussion will be valid even if the sphere is an insulator. We 
assume that the permeabilities of the sphere and the fluid are equal. 


**These conditions imply that the tangential components of electric 
field will be continuous across the sphere and that E = 0 at infinity. 
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and obtain a set of equations from (ll), (12), (15) and (16) with 
appropriate boundary conditions. It should be noted that some of the 
complications which arise from boundary conditions in magneto-gas- 
dynamics are conveniently absent in this problem. 
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III . Solution of the Resulting Equations 

Zeroth Order Equations 

(i) curl H = 0, div E =0, that is the magnetic field to 
this order is independent of the fluid velocity. The same equations 
hold inside the sphere and the boundary conditions are that H q = i at 
infinity and it is continuous across the surface of the sphere. The 
solution is, therefore. 


H 


— o 


= i 


(19) 


(ii) The electric field E satisfies curl E =0, div E = 0, 
' — o — o — o 

both in the fluid and the sphere. The continuity of the tangential 
electric field on the surface of the sphere and vanishing of E q at 
infinity implies 


E = 0 . (20) 

— o 


(iii) The zeroth order flow fields and pressure fields 
are unaffected by the magnetic field and are given by* 

(V 2 + h 2 )v o = - i(^ 2 ) V s 0 - 6 R e v o x i (21) 


O' it 

*The boundary condition v = u jL=-e ion the surface of the 

s a 

sphere and the linearity of the system of differential equations imply 
that the only time dependence of the physical variables is a factor e 1 ^. 


12 


and 


div v = - i s 

-o o 


(22) 


"where 


k 2 


3 Re 

U + 5 R C 2 
e 


(23) 


h 2 


= - i R . 

e 


(24) 


As in the non-rotating case, we have here axial symmetry in that 

the physical variables are independent of the azimuthal coordinate cp. 

However, the azimuthal components of velocity and magnetic field will 

he non-zero here in contrast to the non-rotating case in which v^ and 

H were zero. 

9 


From (21) and (22) it follows that v = 0(l), v = 0(l) and 

cx or 


v 


ccp 


=0(6) and hence 


(V 2 + k 2 ) s q = 0(6 2 ) . (25) 

It follows that the coupled equations (21) can he simplified in this 
case. In fact, in view of (17) and (25), v^ x and v satisfy identically 


*(r, cp, *) are the cylindrical coordinates. 
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the same system of equations and "boundary conditions as in the non- 
rotating case. The solution is therefore [9, p. 656] 


(i) s q = A x fi (k R) 


<“> 2o 



( 26 ) 


where the first term in (ii) is a particular integral of (21) and (22) 
and the last term 

v c = B[ 2 f (hR) ^ x - h 2 R 5 f (hR) V | ] ( 27 ) 

is the complementary function. (R, 0 , cp)* are the spherical polar 

coordinates. The function f (£ ) , defined as 

n 

• (28) 

represents spherical waves with rapidly diminishing amplitude — the one 
with argument C = h R denotes boundary layer effects (viscous) while 
the one with argument C = k R essentially represents compressibility 
effects (sound waves). The constants A and B are given by 


* x = R cos 9 , r = R sin 0 . 
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- i k 2 hfi(h) u s , 

2kf 1 (h)f 1 ’(k) + hfx'Ui) Lkfi’Ck) + fi(k)J 


kfl'(k) Uo 

3L2kf 1 (h)f 1 '(k) +hfi'(hHkf 1 '(k) + fi(k)jj 


(29) 


Now v^ can be determined in a straight-forward manner from (2l). In 

fact, exactly similar arguments can he applied while determining 

Hi and Vi to show that only azimuthal components and are 

affected hy rotation. In view of this decoupling of the system of equations, 

the rotational effects, which only generate azimuthal components of 

velocity and magnetic field, are discussed in Appendix and henceforth, 

we would only consider non-rotating case. 


First Order Equations 

(i) The first order magnetic field Hi satisfies 


V 2 Hi 


i s 1 


o — 


d Vo 

TT ’ 


div Hi = 0, 


whose solution is 


Hi 





-v V C 
o —O 

3 x 


(30) 
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where the first term is the complementary function 

Hi 0 = J BeL g°° 8 > (30a) 

and the remaining terms represent a particular integral. 

On the other hand, inside the sphere, curl curl Hi = 0, div Hi = 0, 
which shows that the sphere acts like an insulator up to this order. 

It follows that 

Hx = D v [R 2 P 2 (cos 9)] (31) 

inside the sphere. The constants C and D in (30a) and (31 ) axe determined 
by the requirement that Hi is continuous across the surface of the 
sphere . 


(ii) Finally, the first order flow fields and pressure fields 
are given by 


cfc- 


= HtP'Hi* - -gi] - ’ Si 


dx 


ka 


(52) 


div Vx 


- 1 Si 


from which it follows that 
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(v 2 + k 2 ) Sl = 3 k 2 [| 2 ||° + i + 


3 v 


ox 


o 3 x 


r . 


(33) 


We have therefore 


s = s c + x ^£o _ Aift „ + vx- 

Sl Sl 2 X 3x 2 X o™ k^h2 3x 


3k 2 ^ v ox C 


(34) 


where the first term 


si = EiR-P^cose) fi(kR) + E 2 R 3 -P 3 (cose) f 3 (kR) 


(35) 


is the complementary function of (33) and the remaining terms are a 
particular integral; and 


v, = v[-l- jv(x -P°) - V ^£o}_ p V s 

- 1 - 1 2 k 2 l ^ 3 x ' h 2 -k 2 3 x 2 J k 2 o r 


+ i^y i + & + ? v 


(36) 


i3 3 Zo ip k 2 +h 2 c . 

2 X 3 x 2 k 2 -h 2 V ox — 


Here the first term 


*v 

ox 

by (27). 


c 


is the x- component of the complementary function 



given 
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v x C = F 1 [2f o (hR) v x - h 2 R 5 f 2 (hR)V | § ] ^ 


+ F 2 [4f 2 (hRMR 3 P 3 (cose)} - 5h 2 R 9 f 4 (hR)V P3 ^ os6 ) ] 


( 37 ) 


is the complimentary function of (32) and the remaining terms are a 
particular integral. The constants Ex, E 2 , Fx and F 2 are determined 
by the requirement that v x = 0 at R = 1. 

The following points are worth noting in our solution: 

(a) The terms (like f n (hR)) containing the coefficients 

B, Fx and F 2 represent the compressibility effects on viscous waves. 

(b) The terms (like f Q (kR)) containing the coefficients A, 

Ex and E 2 represent viscous effects on sound waves. 

(c) The terms corresponding to the coefficients C and D 
represent the potential part of the magnetic field. 

At this point, we would discuss in some more detail the case when 
the fluid is incompressible. 
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IV . Incompressible Case 


The solution in this case can be obtained from the above in the 


limit as c ->cO(i.e., k -* 0), s 0 and 


c s p. 


It is clear that we will still have terms (a) and (c) and the 
corresponding coefficients will be in the limit: 


2f (hj 
o 


[5 + b a l4^J , 


[4 + fi(h)/f 2 (h)] 


C - Cf 2 (h) - 2f 3 (h)] and D - - 


Evidently waves (sound) represented by (b) will no longer be 
present and the corresponding terms from (b) can be obtained by the 
limiting process (38): 


p o = ° 2 s < 


Pi (cosQ) 


2 . r Di 4. oj6 

c Sl L R2 g % 


+ P|| BRfi(hE) - | I^Jp^cose) -t- 


+ L R4 - P t'f ®h 2 R 3 f 3 (hR) 


+ 3 Aji 5C1 

5 R 2 R 4 J 
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v 5 s v s 0 “■ Ai 1 v £ ‘ ( i;f £9) 


v i : 


3 

2k 2 


t 7 C x l°)-h^ v S 0 + 2^ 


2h 2 ds 


es° i - a v (xf (kR)) 


r *- 

J s 


+ P v S X C -> i V j r Dl ^l ( g s6) + D 2 p 3 (cose ) 


R 4 J r 


+ 3 l 3 c V ^ c ° sQ ) + 4 V £ 


R" 3 


£o 


h dx 2 


(40) 


To Is L53 v (f -1°) - R 2 v -I 0 ] 


Hi: 


A s - i + ^4 ?-|° 


7"o S 1 ' ~~"Z\ 

k ^ o — k ' 


v (R 2 P 2 (cos0)) 


_ r 2 v P 2 (cos9 ) 


R3 J 


vhere 


A i - - i h 2 B f 2 (h) , 


3Bh 2 


Dl ' - lOfTThT [ 3 f i 2 < h ) - f 0 <R)fa(h)] , 


Da = - — $ B j_f Q (h) - 5fi(h) + h 2 • 



Next, we can write down the total drag experienced by the sphere* 


X = 


- TT U 

S 



+ = U + s&.r 2 + JS)] 

VfeEe 5 " -v/feEe 



r ' 1 

[21 ^ + 
L. v. 5 



3PHm 1 
5V2Re J 


( 41 ) 


The terms in the curly brackets represent the non-magnetic effect, 
while the last term in each square bracket represents a hydromagnet ic 
effect. The first square bracket gives a frictional force varying as 
the velocity: the first term is the Stokes drag for the classical 

problem in which the sphere undergoes motion of translation; the second 
term is the oscillatory effect in non-magnetic case; the third term 
is hydromagnetic, non- viscous, oscillatory effect and the last term 
is hydromagnetic viscous effect due to oscillation. The magnetic field 
tends to increase the resistance. 


The second square bracket gives the correction to the inertia of 
the sphere. This amounts to the fraction 


x 

2 


+ 






3&rT) 

10 y 


*There will be no contribution of the Maxwell stress to the order 
we are interested in, as there are no currents on the surface of the sphere 
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of the mass of fluid displaced, instead of ■§ as in the frinct ionless, 
non-magnetic case. The effect of the magnetic field is to decrease 
the apparent mass of the fluid. 

Although the derivations strictly do not hold for this case, we 
may consider the value R 00 while at the same time keeping R « 1. 
We then find 


X 


2 " ij i 8 + £ P R. 


i 

5 


m 


U sJ 




a result which can be checked in a straight-forward manner from the 
corresponding hydromagnetic inviscid problem. The non-magnetic part 
provides the inertia effect and the magnetic part contributes to the 
frictional part of the drag. 

2tt 

When the period — is made infinitely long, the drag reduces to 
the Stokes drag in the classical problem in which the sphere moves 
uniformly in a straight line. This is understandable because as 
^ 0, the uniform magnetic field H° should also approach zero in 

view of 


P 





0 ( 1 ) 


and so the magnetic effect is absent in this limit. 
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V. Discussion 

First we will discuss the conditions under which our analysis 
can he expected to give a close approximation to the true flow. 

We have neglected the displacement current in comparison with 
the conduction current, which is justified if 

« 1 , ( 42 ) 

which follows from the order of magnitudes of the various terms in the 
Maxwell equation (6, i) where L is the velocity of light. This is 
satisfied hy all hydromagnetic flow problems in the laboratory. An 
associated condition, that the excess charge may be neglected, is also 
satisfied if (4-2) holds. 

As in the classical problem, the convective terms have been neglected. 
This is justified if* 

o' « a (43) 

(where a is the amplitude of the oscillation and a is the radius of the 
sphere), provided 


*This condition justifies the omission of convective terms both 
in zeroth order and first order approximations as can be seen by the 
orders of magnitudes of various terms involved in the momentum equation. 
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R e = ^ = 0(1) , (44) 

which follows from ( 15 ), remembering that s -* 0, c 2 s ■* p in the Incompres- 
sible case. For common liquids in hydromagnet ic experiments, v = 0(l0“ 6 ) 
in m.k.s. units and hence we must make 

a 2 X = 0(l0 -s ) . (44a) 

Next, in the zeroth order approximation of the momentum equaLion, 
we neglected the Lorentz force J X B = P curl H X H. This is justified if 

e - °< 1 ) to) 

and 

R m = a 2 X p a « 1 . (46) 

In the case of liquid Sodium and Mercury, for example, po = 0(1) 
and so (1(6) follows in view of (44a). Condition (45) implies that 
magnetic pressure is at most comparable to dynamic pressure, which 
holds for weak field, say 50 gauss or less. In short, our analysis 
in valid for the oscillations of small globules in the presence of 
weak magnetic fields. 

The following table gives sane typical values of these parameters 
for Mercury, Liquid Sodium and Saturated Salt Water at 25° C when 
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L, B° 

= 25 gauss, X 

25. 

. - 



Liquid 

Sodium 

Saturated 

Salt 

Water At 
25° C 

j 

Mercury 

R 

e 

T j 

; 3 

1 1 

; 6 

60 

3 

: 25 

25 

2 

R 

m 

j 

11 X 10- 5 

! 

2 X 10" 10 

, 

*85 x 10 “ 5 


Finally, to study the nature and properties of the solution, it 

■will he convenient to introduce the stream functions 5 and ^ 

(where v = - — -=-? v = — -4^, H = - — H = — -=r^) . For the sake 
x r or* r r ox* x r or’ r r 9x 

of brevity, we write the results here in a shorter form: 




i - e sin 2 0 


R + a 2 (R e ,R) e' 


(l+i)yfee/2 (R-l) 


(47) 


+ e iu sin 2 0^ Sd&J + + asfE^R) e -U + iWS572(R-l>j. 


+ e 11 sin 2 e=0B 2 8[Sfi^l + + a a (R e ,K) e" (1+i )^(R-D ] 


(where the terms under the bar represent the non-magnetic part); 


*The wave-like term exp [- (l+i X/Re/2(R-l) ] comes from the functions 
f (hR) defined by (28). 



- 25 - 


5 = e lt cos 0 sin 2 e[b 1 (R e ) + ^ige) + b 3 (R g ,R) e -(l + i)-\^e/ 2 (R-l)J ( 48 ) 


AsE^OOjai--^ , a 3 - 


h 1 6 R « S > 




m a 


m a 


(49) 


ay-'-yiPR - , b x - y= R - 
4 m a * A 4 m a 


and b 2 


. 5_ R a 

20 m a 


If ? is the corresponding function for the magnetic lines of force 
inside the sphere, 

5 = ~ - R R 3 COS0 sin 2 9 e Xt , ( 50 ) 

s 10 a m 9 y 

■which is current-free and hence up to this approximation, the sphere 
acts like an insulator. 

From (47) and (48) it follows that there are three types of terms 
in the solution: 

(l) Wave -like (corresponding to the appearance of a boundary 
layer), effective in a small distance from the surface of the sphere 
and behaving like exp [- (l+i)/y/R e /2(R-l)] . The nature of boundary 
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layer solution is essentially the same in bydromagnetics as in a 
non-magnetic case. The term exp [ - ( 1+i )v^ e /2 (R- 1 ) ] represents a wave 
of vibrations propagated from the boundary of the sphere with the 
velocity J2/Rq, but with rapidly diminishing amplitude, the falling 

■*2tt 

off within a wave length being in the ratio e , or l/535 . The 

linear magnitude ^2/Rg, important in all problems of oscillations in 

incompressible viscous fluids, indicates the extent to which the effects 

of viscosity penetrate into the fluid. In the case of Mercury its 
1 

value is .018 cms, where P is the period of oscillation in seconds. 

i 

For liquid Sodium, the corresponding value is .079 P 8 • 

(2 ) Irrotational (~g“ , COS ^| in —y , and 

/■z\ t, j. u.- ~i /'sin 2 9 sin 2 0cos 2 0 sin 2 0cos 2 0 . 2 q 

(3) Rotational — > g > , sin 0cos0 J. 

In the non-magnetic case, the motion in the flow field is essentially 
irrotational at large distances, and it is an interesting result that 
the magnetic field generates vorticity. 

In the case of a magnetic field, the term . cos ^ in , $ j_ s current free, 
while the term eos0sin 2 0 is responsible for currents at large distances. 

From (^9) it follows that for large distances, the disturbance 
in the present problem differs only in amplitude and phase from the 
one generated by the oscillation of a sphere in a frictionless fluid 
in the presence of a magnetic field. 
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The diagrams 1, 2 and 3 illustrate the effects of various parameters 
on the flow fields . Figure 1 shows how the drag on the sphere varies as 
it oscillates. It also shows the influence of conductivity and viscosity 
on the drag. Figure 2 shows the effect of Reynolds number R g on 
magnetic lines of force. The viscosity tends to contract the lines of 
force . On the other hand. Figure 3 shows the stretching effect of the 
magnetic field on the stream lines (compare Figure 3 (a) and (b ); 

(c) and (d)). Viscosity also tends to produce similar effects on the 
flow (compare Figure 3 (a) and (c ); (b) and (d)^. 
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Appendix 

As pointed out earlier, the only effect of rotation in this problem 
is that azimuthal components of v and H are affected in that they are 
non-zero here in contrast to the non-rotating case in which and H 
are zero. 


To determine v , we use (2l) which gives 


d 2 , d 2 


« ,13 

+ - — - + — -=r- - — - f tr l v 


Sr 2 r dr 




ocp 


5 R v 
e or 


(Al) 


and the boundary conditions are 


= 0 at R =1 and at infinity. 


(A2) 


The solution is (using (26)) 


v ocp = Bl x r fa (hE) + F?k 2 F If 2 


(A3) 


- — ^ ' 6 B h 2 x r fi (hR) 


■where the first term is the complementary function and the remaining terms 
are a particular integral. The constant Bi can he chosen in such a 
way that the boundary condition v = 0 at R = 1 is satisfied. 
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Next, Hj^ can be determined by 


i + S2 + i S _i]H =- 

Ldx 2 dr 2 r "5r r 2 -! "Hp dx * 


(Alf) 


while inside the sphere 


" & +_&!+ li 

_dx 2 dr 2 r dr 


1 . 
r 2 . 


h-, = o 


LJ tp 


(A5) 


and we require the continuity of across the surface of the sphere 


The solution is 


_ Pi(cosQ) + P 3 (cos0 ) , 6E-i d 2 Sn _ 

% " Bs W — B3 53 — + d^ 


3iSB r 2 d 2 j_ d n 

— JF SF + x SL + 1 


r fi (hR) 


(AS) 


+ fi ^ tx r f 2 (hR)] , 


and inside the sphere 


l l 

Hi^ = Ci R • Pi(cosQ) + C 2 R 3 • P3(cos0) 


(A7) 


The first two terms in (A6 ) are the complementary function and the remaining 
terms are a particular integral. The constants B 2 , B 3 , Ci and C 2 are 
determined by the requirement that is continuous at R = 1. 


- 30 - 


Finally, is determined by 


r a 2 , s 2 , 1 s i , 

L dx 2 Sr 2 r 3r r 2 J "^p 


PR v + 6R v i 
e ocp e r 


and the boundary conditions = 0 at R = 1 and at infinity 


The solution is 


v-^ = GiR 2 * P 2 (cos9 )f 2 (hR) + G2R 4 * P 4 (cos0 )f 4 (hR) -f- 


+ 6E e[- as 


Ah 2 -2k 2 ds 0 . x ^Sn __ 
(h 2 -k 2 ) 2 3r h 2 -k 2 SxSr 


_A d 3 so A dfp(kR) 


(h 2 -k 2 ) 2 Srdx 2 h^-k 2 V 5r 


2 S^fr 


h 2 -k 2 3r3x 


1 ^ 1 X r — 


- f 0 (hR) + h¥f 2 (H!)} + 


- 2h2 ( h2 F 1 - w 2 ) X r f! (hE) -f 


+ ^ F 2 x r f 3 (hR) (-Ax 2 + 3r 2 )J - x r f x (hR) 


where the first two terms are the complementary function and the 
remaining are a particular integral. The constants Gi and G 2 are 
determined by the requirement that = 0 at R = 1. 


(a 8) 

(A9) 


(aio) 
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It is clear from the above solution that there will be no 
contribution of Maxwell and Viscous stresses (produced due to rotation) 
on the drag — a result in agreement with our previous conclusions 
(see Singh, 1964 [l2]). 

Since the wave-like terms f^(£ ) are damped heavily, it follows 
that in non-rotating compressible case, the flow will be essentially 
current free and irrotational at large distances. Our solution above 
shows that the effect of rotation will be to produce vorticity and 
current which again agrees with our conclusions in [12]. 

Concluding, it should be noted that the nature of wave propagati cn 
remains essentially similar in both rotating and non-rotating cases. 


- 32 - 

Acknowledgement s 

The author is highly indebted to Prof. J. M. Burgers for going 
through the manuscript and for various helpful suggestions. This work 
was performed while the author held a NASA Resident Research Associateship 
of the National Academy of Sciences - National Research Council. 



- 33 - 


References 

1. R. Hide, "Hydrodynamics of the earth's core", in Physics and 
Chemistry of the Earth , L. H. Ahrens, K. Rankama, and S. K. Runcorn, 
Editors (Pergamon Press, New York, 1956), Vol. 1. 

2. R. Hide and P. H. Roberts, "The origin of the earth's magnetic 
field", in Physics and Chemistry of the Earth , L. H. Ahrens, 

K. Rankama, and S. K. Runcorn, Editors (Pergamon Press, New York 
i960), Vol. 4; , Rev . Mod . Phys . , 32 (i960), 799* 

3- B. Lehnert, Astrophys . J., 119 (1954), 647; , Astrophys . J., 

121 (1955), 481. 

4. S. Chandrasekhar, Proc. Roy . Soc . ( London) , A225 (1954), 173; 

, Proc. Roy . Soc . ( London ), A257 (1956), 476; 

, Proc . Am. Acad . Arts Sci . U. S., 86 (1957), 372. 

5. S. P. Talwar, Jour . Fluid Mech . , (i960), 58I; 

, NASA x-640-64-39 (1964). 

6. J. N. Tandon and S. P. Talwar, Nuclear Fusion , 3 (1963), 75* 

7- S. D. Nigam and P. D. Nigam, Proc . Roy . Soc . (London), A2J2 (1963), 529- 

8. G. I. Taylor, Proc . Roy . Soc . (London), A102 (1922), 180. 


9. H. Lamb, Hydrodynamics , New York, Dover Publ. (1932). 



- 34 - 


10. G. S. S. Ludford, Rev. Mod . Phys . , 32 (i960), 1000; 

, ZAMM , 15 (1963), 9* 

11. K. Tamada, Phys . Fluids , 5 (1962), 817. 

12. M. P. Singh, MASA X-640-64- 142, (1964). 




DRAG VS. TIME 
FIG. I 




TURBED] 
ERCURYI 
16. 2(a) 





FORCE (PERTURBED) AT .95 






STREAMLINES PATTERN (NON-MAGNETIC CASE) AT 
THE COMPLETION OF AN OSCILLATION 
(MERCURY) 

FIG. 3(a) 




STREAMLINES PATTERN (MAGNETIC CASE: B°= 25 GAUSS) 
AT THE COMPLETION OF AN OSCILLATION 
(MERCURY) 

FIG. 3(b) 



STREAMLINES PATTERN (NON-MAGNETIC CASE) 
AT THE COMPLETION OF AN OSCILLATION 

(LIQUID SODIUM) 

FIG. 3(c) 





